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Abstract 

Following work by Khoury and Weltman, we introduce a scalar field 0, the chameleon, which is 
conformally coupled to matter. That is, matter experiences a metric which is a conformal transform 
(parametrized by (p) of the Einstein metric. The effective potential of the field (p is a. sum of its self- 
interaction term and an exponential term due to the conformal coupling. Under certain conditions on 
the self-interaction and the coupling, this effective potential has a minimum which depends on the local 
matter density, as does its second derivative at the minimum. As a result, the scalar field acquires a 
mass which increases with local matter density. 

The field (j) mediates a fifth force which is suppressed in the laboratory and in interactions between 
large bodies like planets, but which may be detectable between small test masses in space. 

In this pedagogical essay, we derive the equation of motion of (j>, discuss chameleon-field cosmology, 
and examine some simple solutions with a view to experimental detection of the chameleon. 
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1 Introduction 

1.1 Why chameleons? 

One of the biggest surprises in modern cosmology was the discovery that the expansion of the Universe is 
accelerating. Taken together with other theoretical and observational evidence [I], this has been interpreted 
to imply that Einstein's cosmological constant A, or some other mechanism for acceleration, must feature 
prominently in cosmology. 

With the Universe spatially flat or very nearly so, the energy density pA associated with a cosmological 
constant (or more generally, with dark energy) is of order Pcriticai = 3i?Q /8ttG today. But in the cosmological 
constant scenario, p\ is constant in time, which means that just after inflation, pA and pr, the energy density 
of radiation, differed by over 100 orders of magnitude [2J. This fine-tuning, and other conceptual problems 
with the cosmological constant, have been discussed extensively in the literature, for example in [UlS]. 

A more general model for cosmic acceleration is a slowly-rolling scalar field, called quintessence (see [4] 
for a recent review). The slow roll means that the field has negative pressure and therefore acts to accelerate 
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expansion. Quintessence can resolve the fine-tuning problem via tracker potentials, couplings between the 
quintessence field and dark matter, or modifications to general relativity. 

Scalar fields, for example the dilaton field, arise necessarily in string theory and are generally coupled to 
matter. But if such a field couples strongly to matter, we should have detected it by now as a "fifth force" . 
A mechanism must therefore exist to suppress the detectable effect of any coupling to matter. 

Khoury and Weltman [S] have suggested a coupling which gives the scalar field a mass depending on the 
local density of matter. In regions of high density (like the Earth), the field has a large mass; in regions of 
low density (like interstellar space), the field has a small mass. In terrestrial experiments, the large mass of 
the field suppresses its interaction with matter. In observations of the Solar System, the action of the field 
is suppressed by a "thin-shell mechanism" , to be discussed. Since the field is able to hide so well from our 
observations and experiments, it is called a chameleon. 

A specific example of a chameleon field arising from string theory is given in [6] , which shows that under 
certain conditions, the radion field of brane theory can act as a chameleon. 

Chameleon fields are not limited to quintessence. Chameleon behaviour can be exhibited in, for example, 
0"* theory, using the mechanism we describe below but with the sign of the coupling constants Pi reversed 

mm- 

1.2 Organization of this essay 

This article is intended to be a pedagogical introduction to chameleon gravity. For a complementary and 
more formal review, see Mota and Shaw 9 . 

We will begin by introducing the action for gravity, matter, and the chameleon field 0, which we will use 
to determine the equation of motion for (p. We will find that we can write down an equation of motion for (jj 
in terms of a bare quintessence-type potential V and a contribution from the matter fields, which will give 
rise to a local relationship between the mass of the chameleon field and the density of matter. 

With this in mind, we will then briefly detour into a discussion of the potential V, finding that an 
exponential potential allows chameleon behaviour while causing cosmic acceleration, but that a power-law 
potential cannot cause acceleration without violating, by several orders of magnitude, existing experimental 
bounds for laboratory detection. 

Next, we will discuss the cosmological evolution of the chameleon field. The principal result of this section 
will be that (f> keeps up with an attractor solution (/imin, analogous to tracker solutions in quintessence, from 
the very early universe until today. 

Then we will derive the equation for the chameleon force acting on a test mass in a static chameleon field. 
A lengthy derivation will give us an approximate solution for static spherically symmetric chameleon fields 
which can be condensed into a useful Maple script. The important result of this section will be a "thin-shell 
suppression factor" which explains why large objects like the Earth don't give rise to large chameleon forces. 

Finally, we will use our solution to show explicitly both why the chameleon field may have eluded detection 
so far and why it may potentially be observed by satellite experiments. 

1.3 Conventions 

In this essay we will set c — 1 and h = 1 and use the metric signature (— , +, +, +). Wc will denote by Afpi 
the reduced Planck mass: 

Afpi = (SttG)"^/^ w 2.44 X lO^^GeV. 

Since we are dealing with what is essentially a scalar-tensor theory of gravitation, we will be making use 
of the terms "Jordan frame" and "Einstein frame". The Einstein frame uses the bare spacetimc metric g^,^, 
whereas the Jordan frame uses a metric g^^l which is conformally related to 17^^ (the index i may be used to 
enumerate matter species, allowing a different coupling for different species). The importance of the Jordan 

(i) . . . 

frame is that in a scalar-tensor theory, gl^l is the metric which matter experiences. In what follows, we will 
make clear which frame goes with which equation. See the first three chapters of pU l for an introduction to 
scalar-tensor theories and conformal transformations. 

For simplicity of notation, we will often write X in place of X^^^^ when we are working with a fixed index 

i. 
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2 The behaviour of the chameleon field 



2.1 The action and its variation for a chameleon field 

Recall the normalised vacuum Einstcin-Hilbcrt action 

1 r . _m: 



EH 



IQttG 



R 



g- 



pi 



R. 



Introduce a scalar field (j) with potential V{(t)) and action 



So far, this is the standard picture. But now we introduce a set of matter fields '0m'' with action 



'S'm — 

which are coupled to by the definition 



d xCy 



('0 o('0 



(1) 



where Pi are dimensionless coupling constants, in principle one for each matter species. Note that £ is a 
scalar density of weight 1 in the language of differential geometry. 
The total action including gravity, 0, and matter is thufS 



5" = 



d^x. 



d^x. 



Ml 1 p 
-fR--{d<t>f 



M, 



'9 ■ 



pi 



R 




(2) 



Variation with respect to <j) allows us to obtain the equation of motion for 

1 d£„ 



5S 




1 dCra dgj^l ^ 
1 dC^, 213, 



9 dg^;^. A//pi 



9df^M^^^^ 

The manipulations of the kinetic term in this derivation deserve careful explanation. On the first line we 
used the symmetry of the metric to get 5 (— iV^c/iV^) — — V^(/)(5V^0. On the second line we used the 
commutativity of differentiation with variation to get SW^(j> — W^Scj). On the third line we integrated by 
parts by applying the divergence theorem to (V^c/)) (50, which is a vector density of weight 1, and 
assuming that — > or (50 ^ at spacetime infinity or boundaries. 
Requiring (55 = we thus have an equation of motion: 



-9 dgi:i ^Vlpi 



(3) 



But the sum on the right depends on both the matter content and 0. We would like to identify and isolate 
the matter content so that we can write an equation of motion with manifest 0-dependence. 



Although we introduced </> as a roUing scalar field which happens to interact with matter, we could alternately have come 
up with S as the action of a Brans-Dicke theory in which the scalar field happens to have a self-interaction potential V {</>). In 
the case that i takes only one value, we could also have written S in terms of the metric glii, making the role of as a fifth 
force clearer. See [101 pp.38-42;66-8], for more on both these points. 
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2.2 Defining the energy density of matter 

In this subsection, we will fix i and write X in place of . 

2.2.1 Jordan frame 

Assuming that the matter fields do not interact with each other [Sj , each energy-momentum tensor 

2 dC 



'9 dhi^ 

is conserved in the Jordan frame [lOl pp. 38-42]. That is to say, 

V^f'"' = 0, (4) 

where V is the metric connection corresponding to the metric g/^^. 

Assume the matter to be a perfect isentropic fluid with p — Wip. Then 

f^^'g^u = -p + ip= - {I- iw,) p. (5) 

This is just the standard form, as in [11' p. 140], written in the Jordan frame. 

Note that since matter of species i couples universally to the metric g^j^, experiments composed of matter 
of species i make measurements in terms of g^i, [12J; in particular, the energy density which we measure is 
P- 

2.2.2 Einstein frame 

Impose, without loss of generality, a Friedmann-Lemaitre-Robertson- Walker (FLRW) background metric. 
The energy density p in the Einstein frame (the frame corresponding to is that quantity, conformally 
related to /5, which obeys the standard continuity equation p oc a^^*-^"*"™'^. 
We will need the following: 

~g^^^ = e2ft^/*V.g,. = e^ft^/^^p' diag (-1, a\ a^ a^) 

= diag(-e2ft*/^^p.,a^a^a2); 

= diag 

Compute the Christoffel symbols: 



~g^''' = diag ( _e-2/3.0/Mp, ^ ~-2^ ~-2^ ~-2 



f 03 — ^02 — f Ji — 2^^" (5a0,l + ffal.O " ff01,a) — 2^^^ (SlO,! + ffll.O — 901a) 

= -a^^ (0 + 2aafi — 0) = a~^a_o; 
f ^3 = ^'22 = f °i = igO" (25„i,i - .gn,a) = ^5™ (2501,1 - .9ii,o) 
= _ lg-2ft0/A/,, (0 _ 2^^^^) = e-2ft^/^^p>aa,o. 

Then use them to expand the conservation equation (jj]) with index /i = 0, remembering when we differentiate 
that we are holding (j) fixed and simply varying the scale factor a: 

= Vyf^" 

rpOl/ I -pO rpav 1 -pz-' fjiOa 



(-P5°"),0 + Se-^'^'^/^^-'aa ora" + Sa^^a (-p5°") 
g-2ft0/Mp, ^ + 3(1 + u;,) 5-ia,o/5) . 
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Thus, multiplying by e2ft"^A^pia3(i+"''): 

= 53(i+™')p,o + 3 (1 + Wi) a3(i+™')-ia,oP 

In other words, in the Einstein frame, the quantity 

p = e3(H-.)ft^/Mpi^ (6) 

is the energy density of matter, because it obeys the continuity equation p oc a-3(i+"'i), (Qf course, if 
(j) <C Mpi, then p « p.) 

Restoring the i subscripts and using equation ([5]), we have the Einstein- frame energy density for each 
species i: 

1 2 dCr, 

^ -(i-3«-i)ft0/M„i 1 2 dCm ({) 
2.3 The chameleon equation of motion 

Substituting equation ([7]) into equation Q , we obtain an equation of motion in which all the (j) dependence 
is explicit: 

V2</> = 1/4^) + ^ (1 - 3u;,) -^p.e(i- 



g3(l+u.O/3.0/Mpi ^ ^ '^^f „(i) 



J^„.p(l-3™,)ft</./Mpi 



We may thus express the dynamics for in terms of a single effective potential. 
Then the chameleon equation of motion in the Einstein frame is simply 

= i4ff.0(0). (9) 

Recalling our choice of metric signature, this means that we expect (j) to seek out minima of Ves- 
Note that in the case of non-relativistic matter, « 0, so the effective potential simplifies to 

i 

We will later be assuming Pi — I, since we wish these coupling constants to take on "natural" values of order 
unity. If all (3i are equal, then we may also assume we are working with a single matter species. 



2.4 The potential and the effective potential 

We wish to choose a bare potential V {(/)) which can give rise to cosmic acceleration today via the slow-roll 
mechanism, as in quintessence |13) . It would be cumbersome to come up with a mechanism that makes 
act as a cosmological constant only today, so we will assume that it has always been rolling down a potential 
slope in, without loss of generality, the positive direction. Thus V (0) should be a monotonically decreasing 
function of (p. 

To obtain chameleon behaviour, we will see below that it is sufficient to assume that V {(p) is of the 
"runaway" form, in the following sense: 
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1. lim^^o V{<j)) = oo; 

2. V{(f)) is C°°, bounded below, and strictly deceasing; 

3. Vcf,{(l)) is strictly negative and increasing; 

4. 1^00(0) is strictly positive and decreasing. 

We concomitantly confine (j) to positive values. Note that although runaway potentials are not found in 
classical physics, they do arise in supersymmetric models; see [l^ for a discussion of such potentials in the 
context of quintessence. 

A convenient example which is often seen in quintessence models (see |15j . for example) is the inverse 
power-law potential (also called the Ratra-Peebles potential) , 

' — 

where M is a constant with the dimension of mass and n is a positive constant. 
Another example is the exponential potential 



y(<A) = MSxp (^— j 



where again M is a constant mass and n is positive. 

The important difference between these potential functions is in the limit lim^^oo ^(0)- The Ratra- 
Peebles potential goes to 0, whereas the exponential potential goes to M*. Therefore, present-day obser- 
vations of the acceleration of the expansion of the Universe give estimates of the value of the constant M 
which differ by several orders of magnitude between the two potentials, as we will see later. This is brought 
to the reader's attention now to avoid later confusion, since Khoury and Weltman [5j and Brax et al. (TB] 
discuss both models. 

If the coupling constants /?; are positive, then although V{(f)) is monotonically decreasing, the presence 
of matter gives the effective potential ([5]) a minimum at a finite field value (f> = 0mini which depends on pi 
(as in Figure[T]). Note that 1^9,0 {4>min) = 0, ie, 

+ E (1 - 3^^) -^p,e(i-3-')'''^™/^'^-' = 0. (10) 
Mpi 



The mass associated with the field (j> is given by 



i 

When (j) = (/>minj we have = m^^^, where 



The reason for this definition, as we will see later, is that m^i^ is the inverse of the characteristic range of 
the chameleon force in a given medium. 

Equation PH]) (see also Figure [T]) shows that if any density pi increases, then (j)„iin decreases, because 
and e(^~'^"''^''''^/^^p' are increasing functions of ^ (it is reasonable to assume w < i). Then we expect 



to increase, because V^^ is a decreasing function of (j). This is obvious from equation pT|) if <C Mpi, 
which we will later show to be the case for the chameleon model we will consider. 

In the case where all the {1 — 3wi) f3i are equal (for example, if we are considering only one matter 
species), it is also easy to show that mmin increases with each pi, regardless of the magnitude of ^min- If 
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Figure 1: A simple example showing the behaviour of physically relevant quantities, with units and constants 
suppressed. Specifically, we plot pe"^ (dashed line), V{(l)) — (/)~^ (dotted line), and their sum, Vca{(p) (solid 
line), p = 1 in the first plot and p = 3 in the second plot. 



(1 — 3wi) (3i = B for each i, where i? is a positive constant, then equation ([TTjl becomes 

= ((/)min) - T-;— ('/>min) • 

JWpl 

Our conditions on V (0) state that V^,^, is a decreasing function of 4> and is an increasing function of 0. 
(/)min is a decreasing function of each pi, and so m,jii,i is an increasing function of each pi. 

When does (p — 0min? We see from the equation of motion ^ that if V^ff.^ = 0, then W'^cf) = 0; for 
example, this is true for static, locally homogeneous solutions. We will also see later that, for the exponential 
potential at least, <p ~ 0min cosmologically from very early times until today. 
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2.5 The energy scale for the potential 

In this section, wc will use the following WMAP 3-year data: [T7] 

~ 0.237, nr>E ~ 0.763; 
Ho « 73.5kms"iMpc"\ 



We will also use 



/^critical 



3Hi 



Pm Pcritical^m 
PDE = Pcriticalf^DE ~ 3.34 X lO""^^ GcV*. 



4.37 X lO-'^^GeV^ 
1.04 X 



2.5.1 The power-law potential 




Figure 2: The value of log^Q M, plotted as a function of n and P, for the power-law potential. Contours are 
drawn for integer values of logj^Q M . 

Suppose V {(f)) — Af*+"(/)^". On cosmological scales, matter is described by a single species with p ^ Pm 
and w — 0. We will denote by 0oo the value of (/)min corresponding to this density. Assume that in the 
Universe at large, (j) = (poo today; we will justify this assumption in section[3l Equation (jlOp . on substituting 
the relation ([6]), becomes 

^"+1 A f ^ ~ 



A. 

Mpi 



4/30<»/Mpi _ 



For this field to account for the dark energy pde, we must now have pde = V (0oo) = M'^^'^(f)J^ . Substitution 
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gives 



Mpi 



" 4/3 V 

where 14^ is the Lambert (or product-log) function [18]. 

For example, supposing n ~ f3 ^ I and using fim ~ 0.237, we get 

0OO ~ 0.477 Mpi w 1.16 X 10^^ GeV. 

Then pDE ~ 3.34 x lO"''^ GeV* gives 

Af=(pDE0^)'/*'+"^«2.O8keV. 
In Figure [21 we explore the influence of n and /3 on the value of M for the power-law potential. 



2.5.2 The exponential potential 

Now suppose V{(t)) — Af* exp (A/" /</>"). This potential is easy to analyse if we assume 3> Af today, for 
then pDE ~ V{(j)) w Af*, and so 

M « (pde)^^"* ~ 2.40 X 10-3 eV. 

There is an important difl'erence between the exponential potential and the power-law potential. For 
large 4> we have 

A.^expf^UAf^ + ^. 

That is, for (j) <^ M, the exponential potential and the power-law potential differ by the constant M*. 
Because this constant comes to dominate in the case of the exponential potential, this model suffers the 
same fine-tuning problem as does a bare cosmological constant A. By contrast, the power-law potential is 
"tuned" by the interaction with matter, suggesting it to be a more natural choice. 



2.6 The chameleon in hiding 

We are now in a position to estimate the range of the chameleon force and to see why the chameleon may 
be obscured from terrestrial observation; we will do this more rigourously in Section 15.21 once we have laid 
the necessary mathematical groundwork in the intervening sections. 

Below, we will assume a single matter species representing the Earth's atmosphere, with patm = lO^^^ g/cm"^ sa 
4 X 10-21 GeV* and u; = 0. 

If (/) 3> Af , then is independent of which of the two potentials we are working with: 

deb j ~ deb V ^ ) ■ 

Therefore we may assume V (</>) = Af^^";/)-" in this subsection for the sake of mathematical simplicity. 
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Figure 3: The value of logj^g "^min in Earth's atmosphere, plotted as a function of n and /3, for the power-law 
potential. Contours are drawn for integer values of log]^Qm,„in. Note that the n and (3 axes are inverted 
here, relative to Figures [2] and [H 



With this assumption, we can determine m,nin exactly for the Earth's atmosphere. Equation (jTU]) becomes 



«+i^4/30„„„/Mp, ^ nAf4+"Mpi 

/3patm 



min g4/30„i„/A/pi(n+l) 



Mpi (n + 1) 

_ Mpi (n + 1) 



4/3 

Mpi (n + 1) 



4/3 



4/3 



Mpi (n + 1) 



nM^+''Mpi 

/3patm 



l/(«+l) 



Then equation pT|) becomes 



mi, 



1/2 



Pl 



We will show later that the chameleon force obeys a Yukawa potential, 



$ (r) = a- 



with A — Current experimental data give an upper bound A < 1 mm for a "strong" (a ^ 1) Yukawa- 

potential force ^19 . This bound translates to m^in ^ lO^^^GeV. We will now explore this bound in the 
(n, /3) parameter space. 

In the case of the power-law potential, we have an exact solution for function of n and /3, which 

we plot in Figure [31 Unfortunately, for this potential we have m^in ^ 10~^^ GeV for all reasonable values 
of n and /3. 

In the case of the exponential potential, we can find function of n and /3 subject to the 

previously discussed approximation with 4> ^ M. We plot this in Figure IH in which the top-most contour 
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Figure 4: The value of logj^Q m„iin in Earth's atmosphere, plotted as a function of n and /3, for the exponential 
potential. Contours are drawn for integer values of logj^g TOmin- 

is TOmin = 10"^"^ GeV. We see that there exists a region in the n-j3 plane which satisfies the experimental 
bound, although it requires both n and /3 to be greater than unity. 

3 Chameleon cosmology 

In this section we will follow the argument of [16] . 

Let us assume the bare potential to be of the exponential form 




with M = 2x 10~3 eV. 



3.1 The cosmological chameleon equation of motion 

Assume a flat, homogeneous, isotropic Universe with the metric 




Then, assuming also to be homogeneous. 





(12) 



which is the usual result for a spatially homogeneous scalar field. 
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3.2 The Priedmann equation 

Suppose the Universe is composed of 4>, pressure-free matter with density coupled to by a couphng 
constant (3, and radiation with density p^.. 

The Friedmann equation is derived from the Einstein equation G'"' — SttGT'"', which can be found, as 
usual, by variation of the action ^ with respect to the Einstein- frame metric 17^,^. But now, for matter of 
species i conformally coupled to </>, we have 

2 3C 



2 dCm dgaf3 



-g dgt^u 

rpp,l'^6f3i4>/Mpi 



Recalling equation we know T"°goo = —p- Using this with the above equation, we have 



~g4ft</./Mp, 



goo 



^g-3(l+™Oft0/A/pig4/3i0/Mp, 
^g(l-3u)i)/3i0/A/pi_ 



The Friedmann equation is therefore 

1 ■ 

2 

We further make the natural definitions 



SH^M^, ^-<P^ + V i^) + pn.e^'t''^-^ + p,. (13) 



and 



Pcritical — 



o = Puil 

Pcritical 

for this cosmology, given the assumption of a spatial flatness. 
3.3 Cosmological evolution of the chameleon 

We know from equation ([9]) that the chameleon seeks out the minimum of Vcs- But that minimizing value, 
0min, now changes over time, as pm is diluted by the expansion of the Universe. We would like to know if 
this process is adiabatic; that is, does the chameleon field 4){t) keep up with the minimum 0min(t) as the 
Universe expands? 

To answer this question, we need to know both how fast (pmXnii) runs away from (j){t) and how fast (j){t) 
can chase after it. Equation (|10p tells us that the characteristic time for the evolution of 0min(^) is roughly 
that of the evolution of pm, i.e., 



Pn 



Pni 



where is the Hubble time. 

The characteristic response time for the field (j){t) can be found by making the damped simple harmonic 
oscillator approximation to equation 
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where the constant lu is given by 



= Kff.^^ (0min) = "^min> 

with the last equahty coming from equation (jlip . The characteristic frequency of the osciUator is of course 
LU = TOmin. If the osciUator is underdampcd, then its characteristic response time is the inverse of the 
frequency, 'm^^. Note that the condition for underdamping is that 2mijiin > 3H. 

In order for (j){t) to keep up with 0niin {t), we must therefore have w"?^^ ^ H^^; that is, 

rrimin > H, 

which is consistent with the oscillator being undcrdamped. 

Assuming that f3 is of order unity and n > ^, we have that if (p (t) is slowly-rolling, then rrimin ^ H 
from the earliest times (ie. the end of inflation) until today (Appendix E]) . If m„iin S> -ff, then 0,nin (t) is 
slowly-rolling (Appendix [B|) . But m„iin ^ H also gives (j) (t) « (pmin (t) , and therefore 4> {t) is slowly-rolling 
too. This picture is consistent: If </> (t) w (f>min (t) at some initial time, then (f) (t) « 0min (^) from that initial 
time until today. 

However, in the future, dilution of pm and will allow V {(f>) ~ to dominate the energy budget. At 
this stage, the Universe will be expanding exponentially as a de Sitter spacetime, so H will become constant. 
However, will decrease quickly as pm decreases, and eventually we will have m ^ H . After this time, 

(p will no longer be able to keep up with (pmim and as the contribution to Vcft (</') from the matter term 
p^P<t> /M-pi flattens out, the dynamics of (j) will be governed solely by the bare potential V {(j)), as in normal 
quintessence. 




3.4 The chameleon in the early Universe 

It remains to discuss the initial conditions mentioned above. We now sketch this out, referring the reader 
to [in] for the details. We assume that the Universe undergoes an initial period of inflation driven by an 
inflaton field with w ~ —1 which is coupled to (j) in the same manner as matter. Then equation ([5]) gives 

l/eff(0) « exp ( ^) + pvace^'^'^/^^^'. 



where pvac is the energy density of the inflaton field. Since pvac is roughly constant during inflation, the 
potential function Ves{4>) is time-independent, so equation is the equation of a damped oscillator. 

Combining equations (jlOp and (|lip and the Friedmann equation, as we did in section 13.31 (or simply 
mapping /3 i— > 4/3 and f2ni f^vac in equation (j31[) ) gives 



Ti^,;,, „ Mr,] 



'mm 



1 + n + n 



M 



48/32r!vac, 



with 

^ ^vac — • 
Pcritical 

Since p and flvac ^-re of order unity, we certainly have, for the discriminant of equation (|14[) . 

> - ■ 48/320vac > 1, 

{2,Hf 9 

so that the oscillator is underdamped in the harmonic approximation. We therefore have that (p is driven 
exponentially towards (/fmin during inflation: 

|0-0min|<e-3^*/2(xa-3/2. 

Note that for the harmonic oscillator approximation to be valid, we must assume initial conditions such that 
m> H initially. 
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At the end of inflation, the inflaton decays into mostly radiation, plus some matter and small excitations 
of the field (fi. Because 1 — Sw = for radiation, it does not couple to (j) and so 0min increases dramatically 
during reheating. We still have m ^ so (j) is driven to and past <?!)inin- But as (f) increases, m decreases 
until friction takes over, at which point (p gets stuck at some elevated value. 

As the Universe expands and cools during the radiation era, matter species decouple from the heat bath 
one by one. When a species decouples, the perfect fluid approximation T^'^ g^^ w —p for matter becomes 
invalid for a time because the decoupling particles have relativistic velocities. This provides "kicks" to 
equation (jl2p which drive (j) back towards (pmin- When the radiation era ends and the matter era begins, the 
driving term in equation \12^ comes to dominate the friction term, so (j) converges to (pmin and then follows 
it, as described in section [331 



4 Chameleons and matter 
4.1 The chameleon force 

The interaction of the chameleon field with matter is encapsulated by the conformal coupling of equation ([1]); 
so is the interaction of the spacetime geometry with matter. Since matter fields tpm couple to g^ll instead of 
to 5^1/, the worldlines of free test particles (meaning particles experiencing only gravity and the chameleon 
force) of species i are the geodesies of g^^l rather than those of g^j, (see also [12)11 
The geodesic equation for the worldline of a test mass of species i is 



xP + TP xt^x" = 0, 



(15) 



where FJ^^^ are the Christoffel symbols and a dot denotes differentiation with respect to proper time f, both 
in the metric. 
Using 

V^Wpi J 
the Christoffel symbols can be determined as foUowsH 



rc 



2 ^ 
1 



iWpi iWpi 



v 



_ 2ft 
ft 



„2ft0/Mpi 



I ik/pi 



Substituting this into (fT5|) gives 



= iP + TP^.xPx'' + {cj,,^5P + - .g'^^^.^M-) i'^i" 

iVJpl 



iP + TP,,x^'x'' 



A 
M, 

ft 



pi 



,xPxP + c^^.xPx'' - g''P<l>,ag^.xPx'') 



iP + TP.^xPx'' + {2cl,,^xPxP + 5""^,.) . 



pi 



^Prom this it is clear that the chameleon force violates the weak Equivalence Principle only if there exist two matter species 
with differing values of Pi. 

^The derivation of this relationship in the case of a general conformal transformation is given in |101 pp.65— 6]. 
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The second term in the above equation is the famiUar gravitational term, while the term in Pi/Mpx is the 
chameleon force. 

We see that in the non-relativistic limit, a test mass m of species i in a static chameleon field (p experiences 
a force given by 

^ = -#V0, (16) 
m Mpi 

as in [5]. Thus, (f) is the potential for the chameleon force. 
4.2 Static, spherically symmetric solutions 

We wish to find time-independent solutions (l){x) for spherically symmetric matter distributions p{r) of a 
single pressure-free matter specie^ in the weak-field limit. 

Assuming that g^u = r]^^, the requirement of time-independence reduces equation ([9]) to 

= yeff> (0) = V.^4> {<t>) + T^pe'3^/*^p'. 

Switching to spherical polar coordinates and assuming spherical symmetry, this becomes 

(P(t> 2d(j) P 



= Vt(Mry) + 4-p('-)^"*"""'^ 



(17) 



Now suppose that we have a homogeneous spherical mass of radius R and density pc sitting in a back- 
ground matter distribution with density pooj with > Poa- Then 



p{r) 

We can define 



Pc r < R 
Poo r > R. 



rric = m (0c) , 



'P=Pa 

Note that (pc < <l>oo, as discussed in section \TM 

We will now discuss solutions to equation (fT7|) . essentially a refinement of the solution scheme given in 

4.2.1 Three classes of solution 



Inside the sphere, equation (|17p drives cj) towards (pc while outside the sphere, it drives (p towards (poo- In 
order to obtain an approximate solution to (|17p. we will assume that outside the sphere, the driving term 
Ves{(p) can be approximated by a harmonic potential, as we discussed in section [531 Then for r > R, 



(P(P 2d(P 2 

+ --r ='^oo(9~ (Poo) ■ 

The general solution to this differential equation is 

0(r)=A + B + (Poo, (18) 

r r 

for dimcnsionlcss constants A and B. Imposing the condition that (p — > (poo as r —^ oo gives i? = 0, so we 
have 

(p{r)=A + (Poo. (19) 



'Clearly, this is identical to the case of a collection of such matter species all of whose coupling constants f3i are equal. 
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We recognize this as a Yukawa potential. We now see why the characteristic interaction range A = ?7i~?j^ 
from section 12.61 was justified, for the resulting chameleon force (I16p retains the exponential suppression 
factor e~'"°°(''~^'. 

For r < R, we will consider two classes of solution based on two different approximations to P7)) . Firstly, 
define i?c to divide the interval [0,R] into two intervals: [0, i?c] on which (j) ^ 0c, and [Rc,R] on which 
4> ^ (j)c. It may be the case that i?c = or i?c = R, so that the interval [0, R] remains undivided. 

We will solve (|17p in each interval, based on two different approximations: 

Approximation 1: (/> 3> 0c- In this case the harmonic approximation to VcS is not valid. But taking a look 
at Figure [U we notice that for > 0min, the bare potential V decays quickly and the term pe'^'^/^^p' 
comes to dominate. In particular, we have (for both the power-law potential and the exponential 
potential) 

(0) « ^Pce^^-^/^p- « pe, (20) 
A/pi Mpi 

since as we have shown, <C Mp\ in all cases today. Equation (|17p now takes the form 

d^cj) 2d(j) /3 



with the general solution 



/? „ .2 , ^ 



'l'i^) = 7r^Py + - + D^c, (21) 
D-Mpi r 



for dimensionless constants C and D. 
Approximation 2: ~ 0c. As with equation p^ . we can use the harmonic approximation 

Kff,0(0)«m?(0-0c), (22) 
but this time we will write the solution as 

cp(r)^E + F + 0C, (23) 

r r 

where E and F are, of course, dimensionless constantslf] 

By smoothly gluing together these solutions and imposing ^ — > as r — s- to ensure continuity of the 
three-dimensional solution at the origin, we obtain an approximate solution to (|17p . There are three cases: 

Case 1: Rc — R. We will call this the "low-contrast solution". Here we use solution outside the sphere 
and (|23p inside the sphere. In order for (r) itself to have a limit as r ^ 0, we must have E — —Fe~"^''^ 
in ([23]), giving 

0(r) = F- + Sr. 



r 

It is easy to check that ^ as r ^ (using L'Hopital's Rule). Therefore, we have the solution 



F 0c r < R 

r 

A h 0OO r > R. 



We have also the gradient, 



F — r < R 

A 2^- r> R. 



^One may wonder why we don't use hyperbolic functions here for simpHcity. The reason for the form we choose for is to 
ensure that the values of the exponents do not become too large for numerical computation. 
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The continuity conditions 



lim (j) (r) — lim (j) (r) and lim — = lim — 

r^R- r^R+ r^R- dr r^R+ dr 



give two linear equations for A and F: 



1 ~ e-2™c-R. 1 
F + (hr^A- 



R " i? 

F [m^R - 1 + mci?e-2™cJ?. ^ g-2m.fij ^ ^ [-rriooR ~ 1] • 

The solution to these equations is 

A = ~ „ (1 - m,R - e-2"^« - mciie-^'"^^) 



mc + moo + TlcC 2mcfl _ 777,^6 ^nic-R 



(1 + mooi?) 



Case 2: Rc — 0. We will call this the "thick-shell solution", in keeping with [5^. We use solution (fT8|) 
outside the sphere and (I2ip inside the sphere. The continuity condition at the origin requires simply 
C = in equation (pij) . so we have 



^ p^ + Dcbc r<R 



(r) 



g-inoo(i'--R) 



r>R, 



and 



3Mpi 



Pc'' r < R 



dr I —rrioore^'^^^'^^^^ — g-m=o(r-_R) 

A 22 _ r > R. 

The continuity conditions give equations for A and D: 



6Mpi " R " 

-pcR A- 



with solutions 



3Mpr" i?2 



A = - p, 



3Mpi 1 + nioaR ' 

1 1\ /3pci?^ 



1 + fUooR 2 y 3(/)cMpi 



Case 3: < Rc < R- We will call this the "thin-shell solution", in keeping with [5]. 

Now the solution is divided between three regions. As in the case Rc ~ R, we have E — —Fe'™"^" 
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Then (j) and d(j)/dr are given by: 



F- + 0c re (0, i?c) 

r 

(^) = i l^P-''^ + - + D^c re (Re, R) 
oMpi r 

p-m^(i — R) 

r G {R, oo) , 
/" 

^rncve e +m^re +e re{0,Rc) 



/3 



C 



6Mp\ 



A- 



-niao (r — R) 



-mao{r — R) 



re (RcR) 

r e (i?, oo) . 



In this case there are four continuity equations: two at Rc, 



F- 



,1 - e 
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C 



F 



Rc 6A/pi 
m^Rc - 1 + TTicRce-'^"'''^'' + e^^^cft 



R. 

P ^ C 



and two at i?, 



(24) 

(25) 
(26) 



The solution to this system is quite unwieldy but is easy to obtain with Maple. However, we still need 
one more equation to tell us what Rc is! 



4.2.2 Determining Rc 

We wish to define i?c G [0, R\ so that of the two approximations (I^Dl) and (|22p. the harmonic approximation 
(|22| is "better" for r e (0, i?c) and the linear approximation (I20p is "better" for r £ {Rc,R). 

How do we define "better"? As (p increases from (pc, the harmonic approximation {(j) — 0c) increases 
without bound from 0. At some point we have {4> — 4>c) — j^—^Pc- But we know that 

Voff,0 ~V^ + -^pc < -r^Pc 
Mpi Mpi 



because < 0. Hence, while the harmonic (Taylor) approximation is the better approximation for 

Mpi ' 



rn^ {(j) — 0c) < -jg—pc once (0 — 0c) > jj—Pc we should switch to the linear approximation. 



Therefore, Rc may be defined by the following procedure: 

• Try Rc = R and compute the low-contrast solution (r) . If 

ml (0(i?)-0c)<^Pc, 

then the solution is valid and Rc = R. 

• Otherwise, try Rc — Q and compute the thick-shell solution (r) . If 

,2 



m," (0(0) -0c) > —Pc 
Mpi 



(27) 



then the solution is valid and Rc = 0. 
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• Otherwise, Rc is defined by the equation 



ml ((/>(i?c) - 0c) = Ti-Pc 

and (r) is given by the thin-shell solution. Substituting the r limit of the thin-shell solution, 

this becomes 



mJ:F 



1 



Pc 



^ F = 



R'c ~ 

PPcRc 

mlMpi (1 - e-2™ci?c 



(28) 



The solution to the four continuity equations gives F in terms of known quantities and R^ so upon 
substituting that solution into (|28p . we have a single equation which determines Rc in terms of known 
quantities but which must be solved numerically. 

4.2.3 Automating the solution 

The algorithm described above is straightforward to follow and has been encoded as a Maple script (Listing 
[1]). This allows us to calculate the effect of the chameleon field for any spherically symmetric two-zone 
density model. 

The algorithm can also be used to plot the dependence of the solution type on parameters like R and 
Poo, as we have done in Figure [3 
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Figure 5: The regions on which the thick-shell, thin-shell and low-contrast solutions hold, as a function of 
R and poo- The exponential potential is assumed, with n — 1, P = 1 and pc — lOgcm"'^. 



4.2.4 The thin-shell suppression factor 

Of greatest interest to us is the behaviour of the field (p outside the sphere, for it is there that experiments 
to measure the chameleon force must take place. 

Here we will reconnect with the previous literature by showing that the thick-shell and thin-shell solutions 
agree with [5] in the limit niooR ^ 1 and (in the thin-shell case) F = 0, {R — Rc) /R ^ 1- Note that the 
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low-contrast solution, while not present explicitly in Khoury and Weltman's work, is simply the Rc 
limit of the thin-shell solution. 

In all cases, the exterior approximate solution is 

g-moo(r-ii) 

(r) = A + (t)oo. 

r 

The dimensionless constant A tells us about the magnitude of and thus of the chameleon force. 
In the thick-shell case, assuming rriooR <C 1, we have 

(3 i?3 

-Pc 



3A/pi 1 -I- rriaoR 

In the thin-shell case, Khoury and Weltman assume (j) = (j)c for r < Rc, which translates for us to F 
Then the two continuity equations give 



2Mpi 0e 

On substituting into the continuity equation (1261) . we get 

p ^ p Z^^_A 



Next, substituting into (^5]) gives 



GMpi SMpi i? 2Afpi SA/pi SMpi R 



2Afpr^ 2A/pi 

Ppc 

Finally, Taylor-expand equation ([29]) in Rc about R to get 

/3 3 2Mpi 

^ Pc7:R—S (0OO - <Pc) 



' -ttR^Pc ' 



The approximate external solutions may thus be written in the following form (as they appear in [5]): 



47rMpi V3 '^7 V PPcR'^ 
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The difference between the thin-shell external solution and the thick-shell external solution is the "thin- 
shell suppression factor" [S", equation 16] 3AR/R, where 



AR 
~R~ 



M, 



pi I 



l3pcR^ 



Khoury and Weltman use this factor to distinguish between the thick-shell case {AR > R) and the thin-shcU 
case {AR <C R). Our criterion (P7|) for the thick-shell condition can be translated in terms of AR/R as 
follows: 



ml {m ~ 0c) > ^P. 

, PPcR^ 



2M 



pi 



> 



pi 



Mpil 



1 



PpcR^ 



2 m?i?2 



AR 1 



We may, in general, define a "chameleon suppression factor" 



W = -A 


\ ^ 




47rMpi 







-A 



3Afpi 



Clearly we have w 1 in the thick-shell case and < W < 1 in the thin-shell and low-contrast cases. This 
gives us a means of quantifying the "thickness of shell" of a object in a given background density. 



4.2.5 Generalizations 

There are two important generalizations of our solution to equation (|17p . 

The first is to situations in which pc < Poo- This is the case for modelling the chameleon field in a 
vacuum chamber, for example, to show that the chameleon field cannot be detected in laboratory vacuum 
experiments. To adapt our solution to this case, we could change the solution regions so that Rc > R. That 
is, we would use the harmonic approximation inside the sphere, the linear approximation on the interval 
(i?, Rc), and the harmonic approximation again on the interval {Rc, oo). However, a vacuum-chamber model 
has already been provided in [S], which shows that if the radius of the vacuum chamber is i?vac, then the 
mass of the chameleon field inside the chamber is i?7ac- This will allow us in section [5T^ to use our existing 
solution algorithm with poo tuned so that TOvac — R^c- 

The second generalization is to the case of more than two density regions, for example in modelling 
the Earth with its atmosphere. We can adapt our solution to this case by writing the density regions as a 
sequence of intervals Rca) U {Rc,i, Ri+i) and solving a system of equations for i?c,i and for the continuity 
conditions at all the interval boundaries. However, if all adjacent pairs of intermediate density regions satisfy 
the thin-shell condition, the solution for the whole system is simply the piecewise connection of the sequence 
of two-region solutions. 



5 Detecting the chameleon field 

Reference [5] provides a discussion of several experimental models which demonstrates why the chameleon, 
if it exists, is capable of evading experimental detection to date. We will provide examples here, assuming 
for concreteness an exponential potential for (f) with f3 — I and n = 1. 
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Figure 6: The Earth's chameleon field (above) and a close-up of the region [Rc, R+ {R — Rc)] (below). 
5.1 The chameleon field for the Earth 

Let us model the Earth as a sphere of radius 6 x 10^ m and density 10 g cm^"^ surrounded by an interplanetary 
medium of density 10^'^^gcm^^ |21j . Ignoring the atmosphere!^ our algorithm gives a thin-shell solution 
(Figures [H and [7]) with 



A w -10 
W « 10" 
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Thus, deviations from Newtonian gravity due to the chameleon field of the Earth are suppressed by nine 
orders of magnitude by the thin-shell effect. 



5.2 The chameleon field for a small copper ball 

Now consider a copper ball of radius 1cm and density pc = 8.92 gcm^"^, such as those used in small-scale 
gravity experiments like [52]. 



^The same approximate result is obtained in [E] including the atmosphere. 
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Figure 7: Acceleration of test masses due to the Earth's chameleon field. The length scales match those in 
Figure [6] 

In the atmosphere, with density poo — 1-3 x lO^^gcm^^, our algorithm gives a low-contrast solution 
(Figure [8|) with 

A« -10^ 
W « 0.01. 

In space, with density poo — lO^^'' gcm^^^, our algorithm gives a thick-shell solution (Figure[9]) with 

A« -10^ 
W Kl. 

Finally, we consider the chameleon profile for such a ball inside a vacuum chamber of radius i?vac, as in 
|22j . As shown in [5j, the chameleon field has a mass mvac ~ R^ac inside the chamber. Using equations pO]) 
and we can assign an effective density poo (^vac) for the vacuum chamber and then apply our solution 
algorithm to find the chameleon profile for the copper ball within the chamber. A convenient example is 
Poo = lO^^gcm^'', corresponding to mvac ~ 1-6 x lO^^^GeV and i?vac ~ 1.2m. Then our algorithm yields 
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Figure 8: The chameleon field of a copper ball in air. 
a thin-shell solution (Figure [TO]) with 

A w -10^ 
W « 0.1. 

The coupling strength of the Yukawa-potential interaction between two identical masses is (cf. [5]) 

We see that for W < 0.1, the experimental constraints of Hoskins et al. [551 Figure 13] are satisfied for f3 <1. 
Thus, the chameleon field can be hidden from their experiment. 

However, if the experiment were repeated in orbit with W ~ I, the chameleon field of the copper ball 
would introduce a correction of order unity to its gravitational field: The Newtonian gravitational acceleration 
of a test mass 5 cm from the centre of the copper ball is 

G^^«G-^«1.0xlO-^^, 
(2 cm) 4cm2 s^ 

and the chameleon acceleration of such a test mass is 2.0 x 10^^ ms^^. See Figure [TT] for a comparison of 
these forces at short range. 



Phi (GeV) 




5.3 Satellite detection of the chameleon field 

The SEE Project [23l [24] proposed to test gravity in space, including searching for a fifth force between two 
test bodies. Supposing that the SEE satellite is a sphere of radius 2 m and density 6gcm~'^, it satisfies the 
thick-shell condition for its chameleon field0 Thus, (j) « (j^oo inside the satellite|f| so an experiment taking 
place within the satellite is effectively taking place in space. As we saw above, test masses in space should 
exhibit a mutual attraction due to the chameleon force of the same order as their gravitational attraction, 
assuming /3 ~ 1. The SEE Project, or another space experiment like it, such as APSIS t25j, would thus 
provide a test for the existence of a chameleon field. 

The satellite experiments STEP, GG, and MICROSCOPE plan to test the Eotvos parameter 77 to very 
high accuracy. This provides a bound on violations of the equivalence principle by quantifying the variation 

'^Of course, satellites have thick shells for a large range of radii and densities. 

* Actually, since the experiment will take place in orbit at an altitude of about 1000 km, if>oc will itself be damped by the 
presence of the Earth, but the qualitative behaviour of <j> (r) near the test mass will be unaffected. 
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Figure 9: The chameleon field of a copper ball in space. 

in free-fall acceleration between bodies of differing composition If the chameleon coupling factor /3 

is different for different compositions of matter, these differences will manifest themselves in equivalence 
principle violations which may be detectable by these experiments. 

6 Coda 

6.1 Summary 

We have given a pedagogical overview of the chameleon field. We derived its equation of motion from an 
action and discussed the definition of the matter density and the choice of the potential V, showing that 
a chameleon field with an exponential potential function can account for cosmic acceleration today while 
remaining hidden from experiments. Then we discussed the cosmological history of the chameleon field, 
and finally its interaction with matter today. With the help of a Maple script, we quickly classified several 
examples of the chameleon field due to spherical matter distributions and used them to explain not just why 
the chameleon field has remained hidden from experiments so far but also why it has the potential to be 
detected. 

There are several topics which we did not treat in much detail: String-theoretical motivations for 
chameleon fields; the chameleon's transition from the inflation era to the radiation era; equivalence principle 
violations; and exact constraints on the chameleon model from current experimental data. Many of these 
subjects are discussed in [5] and [T6] . 

SEE or APSIS will provide strong evidence either for or against the presence of a chameleon field. If it 
does exist, many questions will remain: What is the potential function V7 How does the chameleon field 
behave in the presence of strong gravitational fields? 

The chameleon field is an intriguing model because of its promise of explaining cosmic acceleration and 
because of the possibility that it may soon be tested by satellite experiments. Even if it is not detected, 
it may stimulate other ideas for testable alternate theories of gravity, which is an important part of the 
investigation into the nature of dark energy. 

6.2 Outlook 

The chameleon field has the unusual physical property of density dependence. Matter density appears in 
the perfect fiuid approximation of other theories, such as Newtonian gravity and bare general relativity, but 
these theories may be rewritten in terms of enclosed mass or the mass of point particles. In chameleon 
gravity, however, there is no obvious way to escape the perfect fiuid approximation. This must be regarded 
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Figure 10: The ehameleon field of a copper ball in a vacuum chamber. 

as a shortcoming of the model described here, since in the real Universe the notion of local matter density p 
is not well-defined. It is essential that the point-particle regime of the chameleon model be explored to see 
if it is consistent with the perfect fluid approximation. 

Even if we accept the perfect fluid approximation, we must remember that the real Universe is inhomo- 
geneous on scales ranging from the sub-atomic all the way up to the cosmological. Given the non-linear 
density dependence of the chameleon interaction, we must exercise caution in using the chameleon model to 
make predictions for real experiments. A more general solution method or a successful numerical approach 
to solving the chameleon differential equation would allow us to explore the chameleon on a variety of scales 
and to make predictions in realistic scenarios rather than just toy models. 

In discussing the choice of the potential V {(f>), we determined that a chameleon field with the power-law 
potential is ruled out by existing experiments. But the exponential potential we settled on has the same 
fine-tuning problem as a bare cosmological constant A, and so it introduces more problems than it solves. A 
useful chameleon theory should give ^ pm today without resorting to fine-tuning. 

The conformal coupling ^ which we used is not the only possible choice. Different coupling functions 
give rise to different chameleon field theories, in which the thin-shell condition may be easier, harder, or 
impossible, to satisfy [6]. In any case, the behaviour of the chameleon force depends strongly on the choice 
of coupling function, and we have explored only one specific example here. A different choice may alleviate 
the problems mentioned above with the potential V. 

In addition to proposed satellite experiments, chameleon- type modified gravity may have astrophysical 
effects, for example in the gravitational collapse of gas clouds or in the dynamics of planetary ring systems. 
Again, this demands a more general solution method than the approximation we derived in this work. 

The recent review by Mota and Shaw 9 , while not completely answering the above questions, points the 
way forward, exploring a more general class of potentials than we have done here and also constructing an 
"effective macroscopic theory" of the chameleon field for composite bodies made of small finite-size particles. 
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Figure 11: The acceleration of a test mass near a copper ball in space due to chameleon attraction (solid 
line) and due to Newtonian gravity (dashed line). 



Appendix A Proof that mmin ^ H 



We assume that /3 is of order unity, jt. ^ ^ , and (t) is slowly-rolling. 

Some massaging of equation (|lip will allow us to evaluate fn^in/^^- Start by writing out equation ([TOl) 
for this scenario: 

Mpi 



Expand the second derivative of the bare potential 

d 



Vm = ^ (-nM"0-"-iy) 

71 i—n—2 

= [(n + l) + nAf"0 
Then equation (jlip may be expanded as follows 



(n + 1) Af _^ (nAf"0-"-i)^ y 



(30) 



[(n + 1) + nA/>-,;]J nM-4>-J-'V (0^in) + |^p,„e''^»-/^^p' 



1 -\- n + n 



M 



'/'minMpi 



1 + n + n 



'/'minMpi 



/30,„i„/Mp 



AT, 



pi 



A/, 



pi 



We can put this together with the Friedmann equation p3p to compute 



-ff^ Pcritical ^minAfpl 

Mpi 



= 3/30 



1 + n + n 



1 + ?l + 71 

M 



M 



Af, 



pi 



(31) 
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Certainly the above equation implies that 



vmin 

We will consider two cases, 0niin ^ M and (pmin ^ M. 

Case 1: (/)min ^ We have rim > 10~^* from the Planck era onwards [TH], so equation 



(32) 



gives 



m 



3n/317m • 10^° > 3n/3 • 10^ > 1. 



Case 2: 



Af. Now we have V ((/)min) ~ M^, so equation ([50)1 gives 



-„-i 

mill 



mm 



pi 



/3 



nM4+"A.f, 



Then equation (|32p becomes 



iJ2 



> 3n/3r2,„Mpi 



pi 



iAf4+"M, 



Pl 



3n/3 



(3 [M, 
M 



pl 



M4 



n+2 



From the Friedmann equation (|13p . we see that as long as is slowly- rolling, /Ocriticai is a decreasing 
function of time. Today we have /Ocriticai ~ Af**, so from the earliest times until today, we must have 
Pcriticai > M*, and so 



> 3n/3 



/3 ^Mpi 



n+2 



Finally, we need to put a bound on fl^. In the radiation era and the matter era, is an increasing 
function of time, so we'll try to approximate it at the earliest time, Im, found in the (fimm S> M case. 
(pnun also increases with time as the matter density of the Universe is diluted, so let Im be the time at 



which 



M . Equation pop now gives the value of p,n at Im- 



^V^(0min) = TT-Pme' 



ne 



A/pi 



M-'Mpi - 10"^^ GeV 



Recall that today we have pm ~ Pm ~ 10 '^^GeV'*. Assuming that baryonic (and dark) matter had 
decoupled from the heat bath at tM, the relation pm oc a^^ then tells us that the ratio of a today to a 
at tM is roughly 10^°. Using T oc a'^ [20] with T « 2.73 K today, we see T ~ 4MeV at tM, validating 
the assumption of decoupling. All that remains is to compute pr = ^gT^ ~ 10~^° GeV^ to see that 
this era was dominated by radiation, so that ilm w ^ ^ 10~^ at tM- 

We conclude that 



> 3nP 



> 3n/3 



M 



n+2 



-^Q(24n-12)/(n+l)^ 



which is greater than unity for all n > i. 
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Appendix B Proof that 0min {t) is slowly-rolling 



We assume that m^i^ ^ H. 

We must show that <j)f-^^i^^/2 <C V. Take the time derivative of equation ([TO]) : 

Mpi 

Pm jWpl 

=^ <Pmin — -J, / , N 

VcS,<l)<l) (Pmin) 

4in 9il2F2(0^i„) 



Then 



2V^ (0min) ((/)niin) K,ff,00 ('^min) 

9 ^2 (^^.^) 



< 



2 Voff,00 (l/)!!!!!!) V (0min) (0min) 
9 (^^^.^) 



2 "^inin ^ (0min) 1<00 (0min) ' 

where the inequahty holds because V^ff,^^ > V(j,(j,- We can also compute 



1/1<00 y (1 + n + nA/"0-") nM"(/>-"-2y 
n M" 



1 + n + nM^cj)-"^ 0" 
< , — = 1- 



And hence if <C m"^^^, we have 



6^. 
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Listing 1: A Maple script to solve equation pT|) 

Digits := 30: 
with ( linalg ) : 
c := 299792458; # in m/s 

Mpl := evalf(1.2el9 / sqrt(8*Pi)); # in GeV 

M 2e-12; # in GeV 

V := M~4*exp(M"n/phi"n) ; 

Vphi := diff (V,phi) ; 

Vphiphi := diff (V, phi , phi) ; 

gcm3 := 4.31013337073e-18; 

metre := 5.06772886971el5 ; 

printf (" Enter a name for this profile: "); 
profilename :— readline ( terminal ) ; 

printf (" Enter the value for the coupling constant beta: ") 

beta := parse ( readline ( terminal )) ; 

printf (" Enter the value for the exponent n: "); 

n := parse ( readline ( terminal )) ; 

printf (" Enter the radius R of the sphere in metres: "); 
R_metre := parse ( readline ( terminal )) ; 
R := R_metre * metre; 

printf (" Enter the density of the sphere in g/cm"3: "); 
rho_c_gcm3 := parse ( readline ( terminal )) ; 



V. Sahni and Y. Shtanov, | gr-qc/06060"63 | 

C. M. Will, Theory and experiment in gravitational physics. University Press, Cambridge (1981) 
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rlio_c := gcmS * rho_c_gcm3 : 

printf (" Enter the density outside the sphere in g/cm"3: "); 
rho_inf_gcm3 := parse ( readline ( terminal )) ; 
rho_inf := gcm3 * rho_inf_gcm3 ; 

phi_c := fsolve(Vphi + ( beta/Mpl) * rho_c = 0, phi, O..Mpl); 
phi.inf := fsolve(Vphi + ( beta/Mpl) * rho.inf = 0, phi, O..Mpl); 
m_c := aval ( sqrt ( Vphiphi + ( beta ' 2/MpP 2) * rho_c ) , phi = phi_c); 
m_inf := eval ( sqrt ( Vphiphi + ( beta ' 2/MpK 2) * r ho_inf ) , phi = phi_inf); 

lowcontrastfraction := ( phi_inf — phi_c ) / (m_c+m_inf+m_c*exp( — 2*m_c*R)— m_inf *exp 
(-2*m_c*R) ) : 

if m_c'2*( lowcontrastfraction *(l + m_inf*R) — exp( — 2*m_c*R) ) /R)<beta * rho_c /Mpl 

then 

printf (" Low— contrast condition detected . \ n" ) ; 

A := lowcontrastfraction *(l — m_c*R—exp( — 2*m_c*R)— m_c*R*exp( — 2*m_c*R) ) ; 

F := lowcontrastfraction *(l + ni_inf*R) ; 

phi := piecewise (r<R,F*(exp(m_c*(r— R) )— exp(— in_c*(r+R) ) )/r+phi_c ,r>R,A*exp(— 

m_inf *(r— R) ) /r+phi_inf) ; 
closeuprange R_mctrc/20; 

elif m_c^2*(( phi.inf - (l/(l + m_inf*R)+l/2)*beta*rho_c*R'2/(3*Mpl) )-phi_c )>beta * 

rho_c/Mpl then 
printf (" Thick— shell condition detected . \ n" ) ; 
A := -beta*rho_c*R'3/(3*Mpl*(l + ni_inf*R) ) ; 

phi := piecewise (r<R, beta*rho_c*r~ 2/(6*Mpl)+(phi_inf — (l/(l + m_inf *R) + l/2)*beta* 

rho_c*R'2/(3*Mpl) ) , r>R, A*exp(- m_inf * ( r-R) )/r+phi_inf ) ; 
closeuprange := R_metre/20; 
else 

printf (" Thin— shell condition detected . \ n" ) ; 

continuitymatrix := array( [[0 , - 1 , -X,l -Z] , [0 , 1 , , Y-l+Y*ZrfZ] , [ - 1 , 1 , S , ] , [T 

,-1,0,0]] ); 
solutionmatrix := inverse ( continuitymatrix ) ; 
X := R.O*phi.c ; 
Y := m_c*R_0; 
Z := cxp(-2*Y) ; 
S := R* phi_c ; 
T := l+m_inf*R; 

continuityrhs := array( [ beta * rho_c *R_0 ' 3/ (6 *Mpl) — phi_c *R_0 , beta * rho_c *R_0 

^3/(3*Mpl), R*phi_inf-beta*rho_c*R'3/(6*Mpl) , -beta * rho_c *R' 3/ (3*Mpl) ] ); 
ACDF := multiply ( solutionmatrix , continuityrhs ) ; 

R_0_cstimatc := fsolve (ACDF[4] = (beta*rho.c*R_0) /(m_c^2*Mpl*(l-exp(-2*m_c*R_0 

))) , R.O, 0..R) ; 
R_0 := R_0_estimate ; 
printf ("R.O = %e\n" , R.O) ; 
A := ACDF[1] ; 

phi := piecewise (r<R_0 , ACDF [4] * ( exp (m_c*(r-R_0) )-exp(-m_c*(r+R_0) ) ) /r+phi_c , 
R_0<r and r<R, beta * rho_c * r ^ 2/ (6*Mpl)+ACDF[ 2] / r+ACDF[3] * phi_c , r>R,ACDF[l] * 
exp(— m_inf * ( r— R) ) / r+phi_inf ) ; 

closeuprange := 2* (R_metre— R_0/metre ) ; 

end i f : 

printf ("A = %e\n" , A) ; 
W := -A*3*Mpl/(beta*R^3*rho_c) ; 

r := r_mctre * metre : 

plotsetup (ps , plotoutput=cat (profilename , ' phil . eps ^ ) ) ; 
plot (phi , r_metre=0..2*R_metre , title=~Phi (GeV) ~ ) ; 
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plot set up ( ps , plot out put=cat (profilcnamc , ' plii2 . cps " ) ) ; 

plot (phi , r_metre=R_metre— closeuprange . . R_metre+closeuprange , title='Phi (GeV) ' ) 
plot set up ( ps , plot out put=cat (profilcnamc , ' acccll .cps ^)) ; 

plot ( ( beta/Mpl) * d if f (phi , r .metre )*c' 2 , r_metre = . . 2* R_metre , title='Acceleration 
(m/s^2)'); 

plotsctup (ps , plot out put=cat (profilcnamc , ' acccl2 . cps " ) ) ; 

plot ( ( beta /Mpl) * d if f (phi , r_metre ) *c " 2 , r_metre=R_metre— closeuprange . . R_metre+ 
closeuprange , title Acceleration (m/ s "2) ' ) ; 
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